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Abstract 

We study adsorption of liquid at a one-dimensional substrate com- 
posed of a single chemical inhomogeneity of width 2L placed on an 
otherwise homogeneous, planar, solid surface. The excess point free 
energy r]{L, T) associated with the adsorbed layer's inhomogeneity in- 
duced by the substrate's chemical structure is calculated within exact 
continuum transfer-matrix approach. It is shown that the way r]{L,T) 
varies with L depends sensitively on the temperature regime. It ex- 
hibits logarithmic divergence as a function of L in the limit L — > oo 
for temperatures such that the chemical inhomogeneity is completely 
wetted by the liquid. In the opposite case r]{L, T) converges for large L 
to 2r/o, where r^o is the corresponding point tension, and the dominant 
L-dependent correction to 2r/o decays exponentially. The interaction 
between the liquid layer inhomogeneities at — L and L for the two 
temperature regimes is discussed and compared to earlier mean-field 
theory predictions. 

PACS numbers: 68.15.-Fe; 68.08.Bc 
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1 Introduction 



Adsorption phenomena taking place at substrates equipped with geometrical 
and chemical structures have been attracting considerable interest in recent 
years [1-17]. It is stimulated by the development of experimental techniques 
allowing one to imprint substrates with patterns down to nanometer scale 
(see e.g. [3-6]), and a number of theoretical results predicting interesting ad- 
sorption behaviour for systems with broken translational invariance along the 
substrate. These include filling transitions [7-10], discontinuous changes of 
equilibrium droplet shapes as function of its volume [11-12], or condensation- 
type phase transitions [13]. In case of convex substrates, like spheres or 
cylinders, the wetting transitions are altogether excluded by the substrate 
geometry [14-17]. 

The basic theoretical tool to describe adsorption phenomena at micro- 
scopic level is the well-developed density functional theory |TH]. However, 
in numerous applications it suffices to apply a reduced, coarse-grained de- 
scription based on effective Hamiltonian models. This approach captures the 
system's essential features close to surface phase transitions and often allows 
insight into fluctuation effects, which play an important role in the case of 
systems with short-ranged intermolecular interactions (e.g. in the case of 3d 
wetting). As was demonstrated in Refs [19-21], thermal fluctuations often 
modify the critical singularities of line tension associated to three-phase con- 
tact lines (e.g. affect the values of the corresponding critical indices), and 
strongly influence the properties of the contact line itself [21], especially in 
two-dimensional systems. 

Adsorption at planar substrates equipped with a single, well localized 
chemical inhomogeneity has been studied both from the point of view of the 
adsorption layer morphology and the properties of linear excess free energy 
[22-26]. However, most of these studies referred to mean-field (MFT) approx- 
imation, within which fluctuations of the interfacial shapes are neglected. An 
exception is Ref.[2l| where the droplet height and height-height interfacial 
correlations were studied in the asymptotic regime close to the inhomogeneity 
center, far away from the boundaries. The motivation for theoretical inves- 
tigations of this kind of systems is provided by experimental works (see e.g. 
Refs [3-5]), in which solid substrates are imprinted with chemical patterns 
and brought into contact with fluid. Structures of the adsorbed liquid-like 
layers are then investigated as function of thermodynamic parameters and 
the pattern's characteristics. In recent experiments the chemical patterns 
extensions have approached nanometer scale, at which effective interactions 
between the fluid heterogeneities induced by the chemical structure might 
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become pronounced. An expression describing the aforementioned effective 
interactions has been derived in Ref.[2n] within mean-field approach; it will 
be briefly recalled in Section 4. The influence of interfacial fluctuations on 
these interactions and the associated linear contributions to the system's free 
energy has not been discussed so far. 

In this work we aim at filling this gap by considering adsorption in a 
two-dimensional system whose partition function can be evaluated exactly. 
By decomposing the system's free energy into bulk, line, and point terms, 
one obtains the excess point contribution, which is analyzed as function of 
the system parameters. We show, that the system's behaviour is fluctuation- 
dominated and that the point free energy, unlike the one previously obtained 
within MFT, diverges in the limit of large chemical structure's width (2L) 
provided the chemical impurity is wetted by the liquid. The derived effective 
interaction potential is shown highly universal - for asymptotically large L it 
does not depend on any of the system's parameters. 



The system under study consists of a two-dimensional fluid in a thermo- 
dynamic state infinitesimally away from its bulk liquid-vapour coexistence. 
The fluid occupies the half-space above a chemically inhomogeneous one- 
dimensional planar substrate composed of a chemically different domain of 
width 2L placed on an otherwise homogeneous substrate. The domain is re- 
ferred to as type 1 substrate, while the remaining part of the solid substrate 
is of type 2. The substrate extends along the x-axis from —X to X, where X 
is assumed much larger than L, see Fig.l. One of the phases, say the liquid, 
is preferentially adsorbed at the substrate, while the gas phase remains stable 
in the bulk. We denote the wetting temperatures of (infinite) substrates 1 
and 2 by Twi and Tw2, respectively, and assume Twi < Tw2 while the fluid 
temperature T is taken to be lower than Tw2- 

In this paper we consider the statistical mechanics of the liquid-vapour 
interface based on the SOS Hamiltonian 



where it is assumed that the system's state corresponds to a single-valued 
positive function y = l{x) describing position of the liquid- vapour interface 
over the substrate. Here a is the interfacial stiffness parameter, and V{x, I) 
denotes the interface's potential energy favoring its location near the wall. 



2 Model 




(2.1) 
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Figure 1: A planar, chemically inhomogeneous substrate is exposed to fluid at 
its bulk liquid-vapour coexistence. A layer of a liquid-like phase is adsorbed 
near the substrate. The interface at y = l{x) separates the adsorbed liquid- 
like layer from the gas phase which is stable in the bulk. The symbols /i, 
/2, L, denote the liquid layer height at x = —L, +L, —X and +X, 
respectively. 



In the present model V{x, I) takes - for —X < x < X - the following form 

Vix, I) = e(L - \x\)Vi{l) + Qi\x\ - L)V2{1) . (2.2) 

The corresponding partition function with the interface endpoints flxed at 
(— X, Z_) and (X, /+) is given by the path integral 



Zo{l-,l+, X, L) — 




where the factor is included into the Hamiltonian ^2.1^ . In the case of 
periodic boundary conditions /_ = and upon relaxing the constraint of 
flxed endpoint positions, the partition function can be written in the following 
form 

Z{X,L) = j dl_dhdl2Z2il-,h,X-L)Zi{h,l2,2L)Z2{l2j-,X-L) , (2.4) 

where h = 1{—L) and I2 = 1{L) are positions of the interface at x = —L and 
X = L, respectively. The symbol Zi{yi,y2, X) denotes partition function of 
an interface fluctuating above a homogeneous substrate of type i = 1,2, of 
width A with its endpoints pinned at yi and y2. 
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Upon imposing the "initial condition" Zi{yi, y2, 0) = — ^2), each func- 
tion Zi{yi,y2, A) gains the following spectral representation 



Z,{y,,y2,\) = J2'^niyi)i^^^*iy2)e 



(2.5) 



in terms of orthonormal set of functions solving the Schrodinger equation 

1 



V,(y) + E« ij^:\y) = 0. 



(2.6) 



.2a dy^ 

In the present approach the potentials Vi{y) have the following simple form 



00 for y < 
Viiy) = { -Ui/ksT for 0<y<ai 
for Qi < y , 



(2.7) 



where Ui and Oj are positive constants. The above model was analyzed in 
detail by Burkhardt for the homogeneous substrate. In this case the 
partition function takes in the contact potential limit ^ 0, — » 00, 
Uiof = const, the following form 



Ziiyi,y2, A) 



27rA 

-r,e"^^/2"+"'(^^+^^)erfc 



^(yi + l/2) + rjA 




(2.8) 



where erfc[x] stands for the complementary error function 



erfc[x] 



dte 



(2.9) 



The parameter Xj measures deviation of the system's temperature from the 
wetting temperature corresponding to substrate type i.e. Xj ~ T — Twi (see 
Ref.[22|). Within the present framework it is defined by 



Ti 



1 dtp. 



dy 



(2.10) 



where ifj^^ denotes the ground-state wave function for the z-th substrate. 
The considered range of temperatures corresponds to the situation in which 
only the ground states ip^'' (i = 1,2) may remain bound. In this case 
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'^o\y) = v^— 27^6"^'^ and the ground-state energy Eq'^ — ~2^- the wetting 
temperature T\yi the bound state ip^^ disappears in the continuum of the 
scattering solutions to Eq. ()2.6j) and Eq^ — > 0~ - see |2I|,[2H1- 

We now insert the spectral representations (12 .^jl corresponding to 
Z2{l-,h,X — L) and Z2(/2, X — L) into Eq. ()2.4|) . and take advantage of the 
assumptions that X ^ L and T < Tw2- Under these conditions the energy 

(2) (2) 

Eq of the bound state tpQ is separated from the rest of the spectrum and 
the partition function Z{X, L) takes the form 

Z{X,L) = e-24''(^-^) I dhj dl24''>*{h)Z,ih,l2,2L)^i^\h) , (2.11) 

where terms of the order e^o (^-^) have been neglected. The excess free 
energy rj^L, T) due to chemical inhomogeneity is calculated by subtract- 
ing the line terms 2Eq^^L and 2E^\x — L) from the system's free energy 
F = — \nZ{X, L). The subtracted terms correspond to energies of the homo- 
geneous substrates of type i = 1, 2 of length 2L and 2(X — L), respectively, 
with the planar interfaces at their equilibrium positions. Here E^^ = for 
ri > 0. 

This way one obtains the following formula for ri{L,T). 

r roo roo 

7^{L,T) = -2E'i^L- In/ / dhdhi^f^* {li)Zr{hM,'^L)ij'i\h 

Jo Jo 

(2.12) 

which will be analyzed in the following section. 



3 Excess point free energy 



The integral in Eq. (12. 1211 can be evaluated by substituting the expression 
(2.8) for Zi(/i, /2, 2L), changing the integration variables li = I2 = 
and then integrating by parts. As a result one obtains the following formula 



2^^^^L-ln 



1 



(^[n + Taj 



2T^iT^-T^)L + aiT^ + T, 




r 

e 1 eric 



L 



gf||2 gj-f(;. 




n + T2 y ^||2 [n + T2) 

in which the parameter ^ has been identified as the parallel correlation 
length ^||2 corresponding to interfacial fluctuations above homogeneous sub- 
strate type 2 From Eq. ()3.Hl it follows, that at fixed temperature the 
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quantity ri{L, T) is a non-negative, increasing, and concave function of the 
width 2L. Whether t] is bounded from above depends on the sign of ri. This 
issue is discussed in the following subsection. 



3.1 Large L regime 

In the asymptotic regime of large L the behaviour of the quantity rj given 
by (j3.1|] crucially depends on the sign of ri. This is easily understood by 
noting that the argument of the error function diverges either to plus or 
minus infinity depending on the sign of Xi. Also the value of the ground 
state energy e'^^ is negative and equal — ^ for ti < 0, and equal zero in 
the opposite case. To obtain the asymptotic expression for rj, one expands 
the formula 1)3.1^ for L such that |^ ^ 1 and ^ ^ 1 (the case ti = is 
analyzed separately). This way we obtain the following expressions: 

4rir2 i^n t{t^ V / 



2 V 4cr 

ln( -v^7^S.(^)'^') for ri>0 



ln(^Lj for n=0 

- t2t3 

- ^1 ^2 

{r\-T2Y 

(3.2) 

where = ^ provided ri < 0. The central observation regarding the above 
formula is the change of the character of the L-dependent contribution to 
the quantity rj. In the temperature regime T < Ty/i it converges to a finite 
and positive value 2?7o = ln-^^|^^^, which is a symmetric function of ri 
and T2. Actually rjo can be considered as the point tension due to a single 
inhomogeneity at which two semi-infinite, homogeneous substrates of type 1 
and 2 meet 

The finite L correction is negative and decays as (-^)^/^e~^^/^iii ^ where 
(^l|i sets the length scale controlling the rate of the decay. On the other 
hand, in cases ri = and ri > 0, in which becomes infinite, one finds 
the quantity rj to be logarithmically divergent in the limit L ^ oo. This 
divergence is attributed to fluctuations, as it is absent within the mean-fleld 
approach which will be discussed in Section 4. Let us also note that the 
derivative — ^ has the interpretation of an effective force acting between 
two interface heterogeneities concentrated around x = —L and x = L. This 
force is attractive for all of the considered temperature regimes. It decays as 
(;^)^/^e^^^/'^iii for T < Twi, while in the other cases the decay is of the type 

This force has universal properties - it depends neither on the substrates', 
nor fluid's characteristics. More explicitly, it equals — 2^ for T = Twi and 
— ^ in the case T > Twi (in the units where ksT = 1). 
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3.2 Small L regime 



Within the present framework it is also possible to analyze the excess point 
free energy in the regime of system's parameters corresponding to the inho- 
mogeneity's width 2L small compared to the lengthscales The regime 

of small L in the case ri = is explored separately. The coefficient rj vanishes 
in the limit L — > 0, as expected on physical grounds. The formula (3.1) in 
the case of asymptotically small L, has the following form 

+ (3.3) 

a 

where terms of higher order have been omitted. It follows, that r] decays 
linearly to zero as L ^ 0+ and the corresponding coefficient has the same 
form in different temperature regimes. 



4 Mean field theory 



In this section we compare the conclusions of the previous section with results 
for the quantity 7] obtained within mean-field approximation based on the 
interfacial Hamiltonian 



n[i] = j 



2\dx/ 



(4.1) 



where E is the liquid- vapour surface tension, and uj{x, I) denotes the effective 
potential of interaction between the substrate and the liquid- vapour interface. 
It is approximated by a piecewise constant function of x [231, [2S1 

uj{x, I) = e(L - \x\)uji{l) + e{\x\ - L)uj2{l), (4.2) 

where - in the case of short-ranged intermolecular interactions - Ui (z = 1, 2) 
has the following structure: 

uj,{l) = ne-''^^ + 6e-2'/^^ + ... . (4.3) 

Here denotes the bulk correlation length in the adsorbed liquid phase and 
the positive parameter h takes account of repulsion of the interface from the 
substrate at short distances. For simplicity, h is assumed to have the same 
value for both substrate types. The mean-field analysis of Eq. p.ip . see Refs 
[2Sl,[2n|, leads to the following conclusions: 

- The mean-field expression for rj converges in the limit L — > oo to a finite 
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positive function 2riQ, where the point tension rjo depends on the model pa- 
rameters Ti, r2, S, b. 

- The dominant L-dependent correction to 2?7o is negative. For temperatures 
T < Twi it is of the order e"^^/^!!"^^"^ with amphtude depending on all 
system's parameters; here i\\\MFT = In the remaining temperature 

regimes, i.e. for T^yi < T < T\y2, the corrections are of the order j^, and 
their magnitude is governed by a universal prefactor 

Comparing these predictions to formula (j3.2j] one concludes, that the 
singular behaviour of r] occurring for ti > and L — oo is due to fluctuations. 
It is also worth noting that fluctuations wash out the dependence of rj on 
the interfacial stiffness. Indeed - the energy scale controlling the L- 
dependent terms in r] within the mean-fleld description is replaced by ksT 
(which was put equal to unity in the calculation leading to Eg. (13 .21) ). As 
regards the temperature regime ri < - one observes the presence of the 
prefactor (j;)^^^ in the expression for i], which is absent at mean-fleld level. 
The length controlling the exponential decay of the L-dependent terms in r] 
is in both cases set by the correlation length describing the range of typical 
fluctuations of the interface over homogeneous substrate type 1. 

5 Remarks on mean-field results for ti < 

In this section we introduce a simplifled interfacial model, within which we 
reproduce results qualitatively equivalent to the mean-fleld predictions in 
the case ri < 0. Physically the procedure proposed below amounts to damp- 
ing the long-ranged fluctuations by an external potential W{x,l), but still 
summing all the relevant terms in the partition function. Within mean-fleld 
approach one does something else, namely takes only the dominant term 
from the partition trace. These complementary approaches yield equivalent 
results regarding the dependence of r] on L. 

Although the calculation presented in this section is limited to the tem- 
perature regime T < T^i, the considerations are on general grounds as no 
explicit form of the interfacial potential is applied. We consider a liquid- 
vapour interface fluctuating in presence of an effective potential W{x,l), 
which is assumed to have the piecewise-constant structure: 

W{x, I) = e(L - \x\)Wi{l) + e(|s| - L)W2{1) . (5.1) 

The potentials Wi{l) {i = 1, 2) have single, well localized minima at l^i, whose 
role is to damp the fluctuations and to localize the average interface's position 
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at /^i and Z^2, respectively, for |x| < L and |x| > L . The system's Hamil- 
tonian is given by ()2.H) with V{x,l) substituted by W{x^l). The partition 
function is written in the form of Eq. ()2.4jl and the spectral decompositions 
of Zi, Z2 are inserted into it. The key characteristic of the potentials Wi{l) 
is that the corresponding Shrodinger operators have only discrete spectra. It 
follows, that the partition function takes the following form 

Z = e-<^(^-^)(l + 0(e(<'-<')(^-^))) j dhj dh[t>f{h)ct>f\h) X 

X '<P^\k)<P^^^\k)e-''''^'' + <P'l\k)<P'l^\k^^^^^ , (5.2) 

where the eigenfunctions ^j*'' and the eigenvalues are computed from the 
corresponding Shrodinger equation 



1 d 



2 



W,{y) + Ef <Pf{y) = Q. (5.3) 



-2ady2 

It is then clear that the quantity 77 



7] = lim [- In Z - 2E^^^L - 2eJ^^ (X - L)] (5.4) 
converges to 

2r/o = jdkj dk<pf\k)<pt^\kW{k)<P^o^\k) < 00 (5.5) 
for L — > cxD and that the dominant L-dependent contribution to 77 is given by 
e-^l<^-<'l^| dkj dUfmf*m^\k)<P^^\k) . (5.6) 
The energy gap e']^^ — E^^ > corresponds to inverse of the interfacial 



correlation length [2i|. All the next-to-leading finite L corrections to t] 



are also easily computed from ()5.2jl . The magnitude of the j-th term is of 
the order g-^^l-^f '"^''l for j = 2,3, ... . 

The above calculation shows that the general characteristics of the L- 
dependent contributions to ri{L,T), namely their exponential decay con- 
trolled by do not depend on the explicit form of the short-ranged interfa- 
cial potential. The only assumptions made above are the piecewise-constant 
form of W{x, I) in Eq.(5.1) and discrete spectra of the Schrodinger operators 
corresponding to Wi{l). The later implies that \imi^ooWi{l) = 00, which 
results in damping the long-ranged fluctuations. As a consequence even with 
fluctuations one recovers the mean-fleld theory result rather than the one 
obtained in Section 3 for the potential in Eq.(2.7) which has both discrete 
and continuous spectrum. These two results differ by the factor L~^/^, which 
is due to the continuous spectrum contribution. 
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6 Contact point height 



A remarkable feature of the mean-field interfacial height l{x = L) at the 
contact point where the substrates of type 1 and 2 meet is that it remains 
finite for all L values and all temperatures. On the other hand, the interfacial 
height l{x = 0) at the domain's center diverges logarithmically in the limit 
L oo provided T > Twi [25]. Whether these conclusions remain valid in 
the considered two-dimensional case, where fluctuations play an important 
role, has to be checked. In the following, within the SOS model introduced in 
Section 2, we construct the probability density -Pl(/) of flnding the interface 
at height / for x = L. Once -Pl(0 is known, the average interface position as 
well as higher moments can be expressed in terms of L and T. 
The probability density corresponding to x = ±L is given by 

^40 jvie~m ' ^6-^) 

which for the case T2 < and upon neglecting terms that vanish in the limit 
X — > oo, can be written as 



JdhJdh4'\h)4'^*{l2)Zi{l2,h,2L) 
Inserting the asymptotic form of Zi{l2, h, 2L) 

Zi{l2,li,2L) = -2rie(-ri)e"i^/'^e"i('^+'^) + 
1 /n-\3/2 

+ V^(l) i^i + ^i')i^2 + r,~') + 0{L-'/') (6.3) 
and performing the integration yields the following formula 



:^e(-rOe-i^^/-e-^^ + ^(g)^/^(/ + r^'Kr, ' - r^')]e 



(6.4) 



The average position of the interface at x = ±L is given by 

l + O(e-^/«ii0) for ri<0 



<^>^-< -^^(1 + 0(1/L)) for n>0. ^^-^^ 



and the standard deviation A{L, T) = a/< >l — < ^ >| converges to a 
finite value in the limit L ^ oo in both cases. Interestingly, for ri < one 
obtains 

lim A{L,T) = lim < / >L . (6.6) 

L— >oo L— >oo 
11 



From the above results one concludes that no sign of critical divergences is 
present in the contact regions corresponding to x ~ ±L. Both the height 
< I >L and the fluctuations' magnitude remain flnite in the limit L — > oo 
at all temperatures. It follows that for non-negative ti the behaviour of 
the coefficient rj is strongly influenced by the system's critical fluctuations 
that occur only far away from the contact region, to which the line or point 
tensions are usually attributed. 

7 Summary 

In this paper we were concerned with properties of the excess point free 
energy 77 for adsorption at a one-dimensional substrate equipped with a sin- 
gle stripe-like chemical inhomogeneity. Our conclusions are restricted to 
short-ranged intermolecular interactions and to thermodynamic states at the 
liquid-vapour coexistence. Under these conditions our predictions are as fol- 
lows. 

• At flxed temperature the excess point free energy is a positive, increas- 
ing, concave function of the heterogeneity's width 2L. It converges to a 
finite value in the limit L ^ 00 provided the stripe remains non- wetted, 
i.e. T < T\Yi. The dominant L-dependent correction is negative and 
decays as L'~^/^e~^''"/^iii, where the decay rate is controlled by the in- 
terfacial correlation length corresponding to type 1 substrate (i.e. the 
stripe). 

• For temperatures Twi <T< Tw2 the excess point free energy exhibits 
a logarithmic singularity in the limit of macroscopic domain sizes, i.e. 
L ^ 00. This singularity is attributed to fluctuations and it is absent 
within mean-fleld approach. 

• The derivative — ^ is analogous to an effective solvation force act- 
ing between the adsorbed fluid's heterogeneities at x ~ —L and x ~ 
L. This force is attractive in all temperature regimes. It decays as 
^-3/2g-2L/^lli ^j^g gg^gg rp ^ Twi, while for T > Twi it is propor- 
tional to 1/L, with universal amplitudes (i.e. amplitudes independent 
of the substrates' and fluid's properties). 

• The asymptotic decay of excess point free energy in the regime of 
heterogeneity's width small compared to parallel correlation lengths 
present in the system is linear in L for all temperature regimes. 
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• The results are compared to mean-field predictions. The basic obser- 
vation is occurrence of the fluctuation-induced divergence of rj in the 
limit L — > oo for Twi <T< Tw2- Moreover, incorporating fluctuations 
washes out the dependence of the force's amplitude on the adsorbed 
fluid's properties described by the interfacial tension S and the bulk 
correlation length ^b- This means that the solvation force asymptoti- 
cally does not depend on any of the system's parameters at all. On the 
other hand, both with and without fluctuations one observes crossover 
of the force from exponential to algebraic decay at T = Twi- 

• We proposed a phenomenological way of obtaining the mean-fleld re- 
sults in the temperature regime T < Twi- It amounts to trapping the 
fluctuating interface in a single minimum potential W{x, I). As a result 
the law describing the decay of the dominant L-dependent contribution 
to T] turns out to be the same for many reasonable choices of W{x, I). 

• In order to characterize the system's state in the vicinity of the sub- 
strate inhomogeneities x = iiL we calculated the probability distribu- 
tion function Pl{1) for large L and showed that the average separation 
< I >L oi the interface from the substrate, as well as the fluctuations 
magnitude < P >l — < I >l, remain bounded in the limit L — > oo 
for all temperatures T < Tw2- This means that for T > Twi the 
transition regions x ~ ±L contain no indication of the critical state of 
the interface around x = 0. This result is consistent with mean-fleld 
predictions. 

Although a direct link to experimental studies performed so far (e.g. Refs 
[3-6]) does not seem possible at this stage, we believe that some of our pre- 
dictions are open to experimental veriflcation. For this purpose let us point 
that the results may be referred to a system in which the solid substrate 
is substituted with a non-wettable Langmuir-Blodgett fllm decorated with 
a wettable stripe whose compressibility is orders of magnitude larger than 
in the case of a solid structure. The adsorbed wetting layer will tend to 
shrink the inhomogeneity due to the effective force acting between the linear 
heterogeneities around x ~ —L and x ~ L. This effect should become pro- 
nounced once the temperature Twi is exceeded. Our predictions regarding 
the properties of the force —g^j^ are in principle open to veriflcation via 
direct measurement of the stripe compression. Certainly a three-dimensional 
system is experimentally more accessible and the influence of fluctuations on 
the quantity 77 is not clear in this case. On the other hand both mean-fleld 
theory and the exact calculation in two dimensions indicate appearance of 
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the long-ranged effective force of universal properties. The same kind of be- 
haviour can be expected in a physical system in three dimensions. 
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